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O ; Abstract 

^T) ' We study two-body D PP decays, assuming that each decay process go through the bare 
^ ■ amphtude followed by elastic SU(3) rescattering, where the bare amplitude consists of (i) the color- 
ed ' allowed and color-suppressed factorization amplitudes and (ii) the short-distance weak annihilation 
' amplitudes. We have performed the fit on 14 branching ratios of D ^ PP decays in the 
0^ ■ formalism of the above mentioned model. The final state interactions can be well accounted for 
Q>^ ' by the short-distance annihilation topologies and SU(3) rescatterings. The two SU(3) rescattering 
^ ■ phase differences are 6 = 627 — ^ —46° and a = 627 — 5i ~ —21°, where 627, Sg and 61 are the 
CD ' rescattering phases of final states corresponding to the representations 27, 8 and 1, respectively. We 
^ ■ find that the K^lC decay occurs mainly due to the nonzero short-distance weak annihilation 



effects, originating from SU(3) symmetry-breaking corrections to the distribution amplitudes of the 
final-state kaons, but receives tiny effects from other modes via SU(3) rescattering. Our results are 



ix. 

' in remarkable accordance with the current data. 

^ ■ PACS numbers: 11.30.Hv, 13.25.Ft, 12. 39. St 



1 



I. INTRODUCTION 



It is known that the naive factorization approximation fails to describe the color-suppressed 
D decays. The results can be improved if the Fierz-transformed terms characterized by 1 /Nc are 
discarded jQ]. The short-distance (SD) weak annihilation effects, which may mimic some non- 
resonant final state interactions, have recently been emphasized in two-body B decays Q, l^- 
In D decays, the SD weak annihilation contributions involving gluon emission from the final-state 
quarks, which arise from the {V — A) (V — A) four-quark operators, vanish. Nevertheless, if 
the gluon is emitted from the initial quarks, the SD weak annihilation effects are not zero (see 
the results shown in Sec. IIIB|) and may give sizable corrections to the amplitudes. Such effects 
were first noticed by Li and Yeh 0, 0] and recently discussed in B decays 0, 0] . One therefore 
expects that the SD weak annihilation may play an important role in D decays because the energy 
released to the final-state particles is not as large as that in B decays. Unfortunately, the SD weak 
annihilation topologies, in general, are not calculable in the QCD factorization approach^. 

The color-suppressed B^ D^*'>^{it^ ,r],uj), B^ D^{r]',K^) and B~ — > dec ay in odes 

have recently been observed by the Belle, CLEO and Barbar collaborations |^, 0, [lol. I 111. 11^ . 
These branching ratios (BRs) are much larger than the expectation in the factorization-based 
analysis Using the isospin amplitude analysis of B^ D^it~ , B^ — > D^vr" and B^ D'^tt^, 
one can obtain that the rescattering phase difference of isospin amplitudes ^3/2 and A1/2 is about 
30° It may indicate long-distance (LD) final state interactions (FSIs) are not negligible even 
in B meson decays [l^ . Analogously, larger FSIs could be expected in D meson decays since the 
energy released in D decays is much less than that in B decays as mentioned above. For illustrating 
this point, we perform the isospin decomposition for K^tt~^ ,l{^Tr^ and D'^ K^tt'^ decay 

amplitudes: 



A{D^ - 












A{D+ - 







A 
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Ia 



1/2) 



1/2; 



(1.1) 

where the isospin amplitudes with isospin 3/2 and 1/2 are denoted as ^43/2 and ^1/2, respectively. 
The relative rescattering phase between ^3/2 and A1/2, denoted as (j), satisfies the following relation, 



^ One may introduce the transverse momenta of quarks (k_L) to regulate the end-point divergence, where 
is naturally constrained by the infrared cutoff ~ l/R with R the meson's radius (some other discussions 
can be found in Ref. 0). However, the result may suffer from the gauge problem and is part of higher twist 
contribution. It is interesting to note that in deeply inelastic scattering (DIS) processes, by introducing a 
generalized special propagator |^ for massive quarks , the separation of the hard part (T) from the soft 
part (parton distributions) is manifestly gauge invariant for different orders in 1/Q (twist). An important 
feature of using the special propagator technique is that the contributions should be moved into T, 
such that, after combining with the gluon field A" in T, a covariant derivative of color gauge invariance 
can be achieved and classified as a high-twist contribution. 
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cos(p = ^ ^ . (1.2) 

2V2|^i/2| 1^3/2! 

Substituting the data for BRs of K~ 7r~^ ^ K^tt^ and k'^t:^ modes, which are (3.8 it 

0.09)%, (2.30 ± 0.22)% and (2.82 ± 0.19)% Q, respectively, into Eq. ((ESI), one can obtain the 
rescattering phase (j) ~ 94°, much larger than that in the charmful two-body B decays. The above 
result indicates that FSIs should be significant in D meson decays. 

In this article, we will assume that the FSIs in D ^ PP are described by SD weak annihilation 
topologies and elastic (LD) SU(3) rescatterings. Analogously, the elastic final-state rescattering 
picture has been extended from SU(2)-type to SU(3)-type in B decays [1^ . We presume 
that each D — > PP decay process go through the ^^bare^^ amplitude followed by elastic-SU(3) 
rescattering, where the bare amplitude describing the SD-dominant contributions consists of (i) the 
usual factorization amplitudes of color-allowance and color-suppression, which can be calculated 
using the factorization approach, and (ii) the SD weak annihilation topologies (Vl^-exchange or 
14^-annihilation) which present the endpoint singularities are regulated by introducing the complex 
phenomenological parameter Xa 0| in the QCD factorization approach (see the detailed description 
in Sec.lTTBl). 

Interestingly, the SD weak annihilation amplitudes are dominated by the topologies of gluon 
emission arising from the initial-state quarks of the weak vertex, while the total amplitudes vanish in 
order of if the gluon is emitted from the final-state quarks. On the hand hand, the elastic SU(3) 
rescatterings are mainly generated by gluon exchange between the final-state mesons. Therefore, it 
could be expected that the possible double counting is negligible between the two possible sources 
for FSIs. We will give a detailed discussion for possible rescattering sources in Sec. 

We consider the SU(3) breaking effects in the bare amplitude level, but, for simplicity, do not 
distinguish the breaking influence on the two SU(3) rescattering phases, defined as 6 = 827 — 5s 
and a = 627 — 5i. In other words, in description of decay amplitudes, masses vary according to 
SU(3) breaking, and meson productions differ in strength as reflected in the decay constants and 
form factors. 

The rest of this article is organized as follows. In Sec.^ neglecting elastic SU(3) FSIs, we first 
sketch the factorization amplitudes as well as the SD weak annihilation contributions in two-body 
D decays. Sec. IIIII is devoted to the formulation of SU(3) rescatterings. We give the numerical 
analysis in Sec. IIVI The discussions and summary are presented in Sec. The detailed results for 
the factorization amplitudes, SD weak annihilation amplitudes and tensor approach for the SU(3) 
final-state decomposition are collected in Appendices El El andjUl respectively. 



II. THE BARE AMPLITUDES 

Here we present factorization and SD weak annihilation amplitudes for D — > PP decays. The 
relevant effective Hamiltonian for the charmed meson decays is 

Heff = ^ VugV:g,{ciOi+C202)+H.C., (2.1) 

q,q'=d,s 
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where Gp is the weak couphng constant, and the current-current operators read 

with {uq)y_j^ = u'-/'^{l — 7^)(?. Vuq and V*^, are the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements given by 

/ V^d Vus Vub\ f 1 - f A AX^ip - ir]) \ 

Vcd Vcs \= -X 1 - f AX^ (2.3) 

\Vtd Vts Vtb J \AX^{1- p-i7^) -AX^ 1 / 

in the Wolfenstein parametrization. 

The two ingredients of the ^^bare" amphtude for describing the decay processes are (i) the 
factorization amphtudes, which are made of the color-ahowed external Ty-emission tree amplitude 
(T) and/or the color-suppressed internal Vl^-emission amplitude (C), and (ii) the weak annihilation 
amplitudes which consist of VF-exchange and/or ^^-annihilation topologies. 



A. Factorization amplitudes 



Taking K tt^^K^tt^ as examples, the factorization amplitudes can be written as the 

following general forms: 



-r^"^" = ^VudV:,a,iUml-ml)F^^{ml), (2.4) 
^ ^VudV:,a2ifK{ml-ml)F^-{ml\ (2.5) 

where the superscripts denote the decay modes. Here, the nonfactorizable effects, including the 
radiative corrections to the weak vertex and the spectator interactions, are absorbed into the 
parameters ai^2 which amount to replace (= the number of color) by A^^^ such that 

1 , , 

ai,2 = C2,i + ci,2^. (2.6) 

We have summarized the factorization decay amplitudes in Appendix 1X1 where the physical r]' and 
rj states are related to the SU(3) octet state r/g and singlet state r/o by 



\rf) \ ( cost? — sint?\ / Ir/g) 
jr/') I \ sini? COS1? / \ |r/o) 



(2.7) 



with the mixing angle d = —15.4° (ri |. and 



\rjQ) = ^=\uu + dd + ss) , \ri^) = ^=\uu + dd — 2ss) . (2. 
v3 v6 



Introducing the decay constants /s and /o by 



we have 



/", = ^ sin?? + ^ cost?, = -2Asint? + Acos??, (2.10) 
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(a) ^ (b) (c) (d) 



FIG. 1: Annihilation corrections to D ^ Pi P2, where (a) and (b) correspond to A^, while (c) and 
(d) give rise to A\. 

and 

= ^ cost?- ^ sin t?, = -2Acosi?- ^sini?, (2.11) 

where 

{0\ujf,j,u\v^'Hp)) = if^ioP,^, {0\si^15s\v^'Hp)) = i/^wP/.- (2-12) 
The form factors for B — > ij^'^ transitions are assumed to be 



B. SD weak annihilation amplitudes 



The SD weak annihilation contributions 13, 1^ to D — > P1P2, graphically shown in Fig. ^ are 
represented as 



q,q'=d,s 



(2.14) 



In general, {PiP2\Tb\D) consists of ic/u/pi/pj 61,2, where c contains factors of ibl,ibl/\/2, l/\/6, 
or — 2/\/6, arising from the flavor structures of final state mesons, and 



61,2 = ^Ci,2A{{P2Pl) 



(2.15) 



with the convention adopted here that P2 (Pi) contains a quark (antiquark) arising from the weak 
vertex with longitudinal momentum fraction x (y). Here the basic building blocks for annihilation 
amplitudes originating from operators {qic)v-A{Q2Q3)v-A Stie denoted as A^{^ , where the super- 
script i (/) indicates gluon emission from the initial- (final-) state quarks in the weak vertex, given 
by 



A{iP2Pl] 



TTOs I dxdy <j $P2(x) ^>Pi(y) 
0. 



1 



with r^' being defined as 



y(l - xy) 



2m|. 



+ 



1 

x'^y 



+ r^^r^^1>^P^{x)^Pp^{y) 



xy 



(2.16) 



and ruq^^q^ the current quark masses of the meson constituents in the MS scheme. The relevant two- 
parton hght-cone distribution amphtudes (LCDAs), up to twist-3, of a hght pseudoscalar meson 
P are defined as 0| 

Jo 



{P{p)\q2{z2Wul5qi{zim = ifpfip{p,z,-PuZ^.) ['dxe'^^P-^'+^P-'^^^^, (2.18) 

JO o 

where z = Z2 — zi, /xp = m?p/{mq^ + mq^), fp is the decay constant, and x (or x = 1 — x) is the 
cohinear momentum fraction carried by the quark q2 (or antiquark qi). Here and below we do not 
explicitly show the gauge factors 



Pexp 



igs / dt (zi - Z2)^,A^'{tZl + (1 - t)z2) 







(2.19) 



in between the quark fields. The leading-twist LCDA $p(x) is of twist-2, while $p(x) and <I>p(x) 
are of twist-3. LCDAs appearing in the calculation of weak annihilation contributions are in the 
form of 

{P{p)\q2A^2)qiA^im 

= ^ j\x e^(-P-2+xp..i) 1^^^ _ ^^^^ I^^P _ ^^^^^u j I (2.20) 

Neglecting three-particle contributions, the twist-3 distribution amplitudes in the asymptotic limit 
are related to each other by equations of motion, so that 

^Pp{x) = l, = (x - x)c^?, , iM^ = (xx)cI>^. (2.21) 

Using the above simplification, one can get the corresponding projector of Eq. 1)2. 2U() in the mo- 
mentum space 0, 0, Q| 

= — ( ^75 ^P{x) - ^P75 -^^T^ 



M^p = ^ ( ^75 ^p(x) - fiPTo ^p(^) ) ' (2-22) 

a/3 



and further obtain the basic building blocks for annihilation amplitudes given in Eq. ()2.16)1 . where 
the momenta of the quark qi and anitiquark q2 in a meson are parameterized as 

k>i = xEn''_+k''^ + ^n^, k^ = xEn^_-k^^ + ^n'X, (2.23) 

respectively. For simplicity, we have introduced two light-like vectors = (1, 0, 0, — 1), n!|^ = 
(1, 0, 0, 1). If neglecting the meson mass squared, we have = En^ where E is the energy of the 
meson. We refer the reader to Refs. for the detailed technique of calculating weak annihilation 
contributions. 
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The LCDAs normalized at the scale /i can be expanded in Gegenbauer polynomials of forms 



$p(x,/^) = Gx{l — x) 



n=l 



n=l 



^p{x,iJ,) = 6x(l 



^(/.)C(i/2)(2^_l) 



n=l 



(2.24) 



(2.25) 



(2.26) 



In the numerical analysis, we truncate the expansion of $p at n = 1 and just take the asymptotic 
approximation for $p and ^>p. Note that af is nonzero only for the kaon. For the kaon containing 
an s quark, we have the replacement x x in Eq. (|2.24() . The annihilation corrections to — s- 



as an example, thus read 
7O 



ArJK^K 



.Gp J, J, 



2 Cp 



dxdy 



+VudV*d / dxdy 



Jo 



^^{x)^Ko{y) 
$^o(x) $-^o(y) 



1 



-— + -J— 1 + (r 

xy) x'^y 



K xy 



+ 



^^/D/i^Tra^ciK^VcsSGaf (4Xa + 33 



iV2 



xy) 

-47r2), 



1 

xP'y 



+ r 



$P (y) 



K 



2 



(2.27) 

where use of VudVcd ~ "^sV^, has been made, and dzj z Xa has been used to parameterize 
the logarithmically divergent integrals Q, 0| , which can be regulated by including the transverse 
momentum of the quark in the end point region of integrals, but however may suffer from some 
theoretical problems (see discussions in the introduction). It is interesting to note that Ap^(K^J<^) 
is proportional to . As will be seen in Sec. IIVI the magnitude of has a large impact on 
the K^¥^ branching ratio. Two remarks are in order. First, the simplified form of the 

projector in Eq. ()2.22() cannot be justified if considering higher Gegenbauer moment corrections to 
$p and <I>p. We have checked that the amplitude corrections due to a^'^ and a^''^ are numerically 
negligible if the magnitudes of a^'^ and a^'" are not too large. Second, we do not consider , 
since distinguishing a2,of , and is not numerically significant in the present study and, 
moreover, partial effects due to can be absorbed in Xa- The detailed expressions for SD weak 
annihilation amplitudes are collected in Appendix IbI 



III. SU(3) RESCATTERINGS 

From the isospin amplitude analysis of K~'k~^ , K^tt^ and — > K^7r~^, as discussed in 

Seed we know that the LD FSIs effects may be significant in D meson decays. Considering elastic 
SU(3) rescatterings in D decays, for instance, — > K~tt^, — > K^-ir^ and 'K^Vs can be 
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generated via 



K--K+ 
K V8 



K--K+ 
K V8 



1)0 



K r]s 



K-TT+ 



Taking into account elastic SU(3) FSIs, the decay amplitudes Af^^ are given by 2^, 21, 

aFSI _ cl/2 A bare _ /jtTc1/2 tt\ a bare 



I 



(3.1) 



where S is strong interaction scattering matrix, and A^^^'^{= Ap'^ + A.J'^) are approximated in 
terms of the factorization and SD weak annihilation amplitudes. Note that S is unitary. The 
SU(3) final-state rescatterings for D — > P1P2 are described by the product 8 8. Since the 
P1P2 states obey the Bose symmetry, only the symmetric states given by the representation 
36(= 27 © 8 © 1) in 8 © 8(= 36 © 28) decomposition are relevant, whereas states given by the 
representation 28(= 10 © 10 © 8) vanish. 

In the present study, we will use 527, ^8 and 61 to stand for the respective rescattering phases 
of 27, 8 and 1 states. The detailed derivation for U matrices and the corresponding SU(3) eigen- 
amplitudes is exhibited in AppendixO Thus the S^/^ matrices and decay amplitudes can be recast 



into the following 5 subsets (see also Refs. |l6l. l22l|^: 



subset 1 {K TT"' 



^0 — 

K TT^ 



K ?78 rescatterings), 

/ 2+3e-^^ 3(1- 



cl/2 



■1(527 



5 

3(l-e-'*) 


5V2 


5V2 
x/3(-l+e-'*) 


5^2 


10 


10 




) y3(-l+e-'*) 




5V2 


10 


10 



(3.2) 



bare 



( J^^'- + 
K 7r+ 

/(bare 

^bare 
K ris 



V 

subset 2 {K^TT^ — K^ir^ — K^rjs rescatterings), 

^.1/2 _ ^.1/2 



A bare 



K +TT 

/(bare 



(3.3) 



(3.4) 



(3.5) 



A bare 

subset 3 {K^7r~^ — K~^7r^ — K~^r]s rescatterings), 



■i527 



2+3e- 



3(1- 



v/3(l- 



3(1- 



5V2 



5\/2 
v^(l-e- 



7+3e- 



5\/2 
V3{-l+e- 



5V2 



A bare 



10 

V3(-l+e- 
10 

/ A bare 
/(bare 



10 



9+e- 



10 



(3.6) 



(3.7) 



lare 
K+ris 
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subset 4 (tt'^tt^ — tt'^ijs — K~^K rescatterings) , 

(I 



(TTTr)"* 





Vo - 



3+2e- 



V6(l- 



\ 

v/6(l-e-"') 
5 



* bare 
-^(7r7^)^ 



5 5 
/ /(bare \ 
yj^bare 

TT+TJS 

/ibeire 



(3.8) 



(3.9) 



• subset 5 (7r"'"7r — tt^tt^ — r/8??8 ~ K^K — K^K^ — n^ris rescatterings) 

Sl/2 —1(527 

(tttt) 



5e- 






5e- 


^''+8e-»''-13 


5e- 


""'-8e"''^ 


+3 


5e-»<'-4e-*''-l 


5e-i<^-4e- 


i(5 -[^ 


5e- 


20 


■5-13 


5e- 


20v^ 

»<'+8e-"5+27 


5e- 


20v^ 

-"'-8e-"5 


f3 


20 

5e-»<'-4e-»''-l 


20 

5e-*<^-4e- 


i(5 


5e- 


20V2 

-'<^-8e- 


**+3 


5e 


40 


5e- 


20v^ 

"'+8e-''^^ 


-27 


20V2 

5e-i<T_,_4g-i5_g 


20V2 

5e-»<^+4e- 


»<5_9 




20^2 






20v^ 




40 




20v^ 


20V2 




5e- 


-<<T_4g- 


i<5 


5e- 


-»o-_4g-i«_]^ 


5e- 


-»<'+4e-i'5 


-9 


5e-'<'+8e-*''+7 


5e-»<^_4e- 


i<5 




20 






20\/2 




20\/2 




20 


20 




5e- 


-i<'-4e- 




5e- 




5e- 
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5e-io-_4g-i<S_]^ 


5e-*''+8e- 






20 







20a/2 






2OV2 






20 

4V3(e-<''-l) 
20 


20 

4V3(l-e" 
20 










4x/3(e-'*-l) 
20 

20 

4(2e-'''+3) 
20 



(3.10) 



A bare 
-^(7r7r)0 



(3.11) 



larc 

/I bare 

/I bare 
^K+K- 
/ibare 

^bare 

where b = ^21 — a = ^27 — <5i, and we have included the identical particle factor l/\/2 in the 
amplitudes -A^o^^o and -4j^gjjg. Here S-*^/^ matrices have been factored out an overall phase e"'^^ since 
only phase differences affect physical results. Note that we do not list — A^tt''", which does 
not belong to any above subset, i.e., does not rescatter with other PP modes. Note also that in 
the subset 4, — tt+tt^ does not rescatter with other modes, too. 



IV. RESULTS 

In this section, we will first introduce the relevant parameters in the fit and then give the 
numerical results together with a brief discussion. The 2-body D meson decay rates are given by 

V {D ^ P^P^) = ^\A^^^\\ (4.1) 

where ^ is the center-of-mass momentum of decay particles. In the numerical analysis, we perform 
the best multi-mode fit for measured branching ratios, defined as 

X^ = y,f^)'. (4.2) 



Ax,; 
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where and Xi it Axi denote the theoretical results and measurements, respectively. On the 
theoretical side, input parameters relevant for our numerical analysis are listed in Table [isl . 
3, Q, Esl- As listed in Table HB we take the current data for the 14 Ktt, tttt, KK, Kr)^''^ 
and 7rr/(') BRs as inputs. The modes involving 77 or r]' are related to r/g and rjQ via the mixing 
angle i}. The SU(3) FSl picture is not suitable to be extended to the U(3) scenario since Uyi(l) 
symmetry is broken by anomaly, i.e., r]' is not a Goldstone boson. The weak annihilation effect 
for SU(3) channels is parameterized in terms of Xa, while that for decay modes involving 779 is 
distinguished to be X'^. However we do not distinguish because it is numerically small, 

as seen in our analysis. The scale for the factorization amplitudes is taken to be /i = nic, i.e., 
l/N^^ = \/N^^{mc), while the scale for SD weak annihilation amplitudes is 1 GeV. We use the 
world average value of Fq^{Q) = 0.76 it 0.03 E^]. For the (f' dependence of form factors, we adopt 
the pole dominance assumption: 

i^o(O) 



(4.3) 



1 - q^/ml ' 

with taking rri-t as the mass of the lowest-lying scalar charmed meson in the corresponding chan- 
nel. The above form is consistent with the recent QCD sum rule study for B — > light meson 
transitions H- We assume = 2.3 GeV (or 2.2 GeV) for F^^ (or F^""). The results for 
fitted parameters, which are (i) two FSl phases, 5 and cr, (ii) the form factor Fq'", (iii) SD weak 
annihilation parameter X^ and X^, and (iv) 1/A'^^''^, are cataloged in Table ITTTI Output observables 
are given in Table Ull The errors of outputs correspond to the variation of F^^{0), while the errors 
due to uncertainties of D lifetimes are negligible. 

The nonfactorizable effects are lumped into the effective number of color N^^, of which the 
deviation from Nc measures such effects. could be complex. However it is assumed to be 

real due to its small value: < —1/15 (~ —0.067) in the fit, consistent with the very earlier 



large- A'^c approach for describing hadronic D decays [1|| . It is interesting to note that we obtain the 
weak annihilation parameter \Xa\ = 3.84±0.06 {\X'^\ = 2A5to^l 2.18±0.19) with a large phase 
(—138 ± 3)° (Y— 138 ± 3)° or (130 it 3)°), compared with the similar parameter \Xa\ ~ 4.5 given 
in B decays 5l. l28l. 129 1 . Note that we obtain a twofold solution for X^. As seen in Table ITTl the 
weak annihilation topologies have a large impact on branching ratios. This analysis gives moderate 
rescattering phases 5 ~ —46° and a ~ —21°. 

One can see from Table UTI that the D PP data can be nicely fitted by the present picture. 



A. D+ TT+vrO vs. D^ 



Consider the ratio 



Vcd 



(4.4) 



T{D+ -^ICtt+) 

The data show Ri = 3.46 it 1.17, whereas i?i = 1 in the SU(3) limit. It is interesting to note that 
both — > TT^TT^ and — > K^vr"*" amplitudes are identical during SU(3) rescattering because 
they do not rescatter with other decay modes. Moreover, these two amplitudes have no SD weak 
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TABLE I: Summary of input parameters on the theoretical side of 

the fit. 



Running quark masses [GeV] and the strong 


couplin 


2; constant 




m,(lGeV) 


m„(lGeV) 


mail GeV) 


a,(l GeV) 


1.35 


0.12 


0.004 


0.009 


0.517 




The Wolfenstein parameter and D-meson lifetimes 


[10-15s] 


A 


t{D+) 




0.2196 


1040 ± 7 


410.3 ± 1.5 


Pseudoscalar-meson decay constants [MeV] 


u 


fx 






fo 


131 


160 


168 


157 


220 ± 20 




The form factor (at q'^ 


= 0) and 7] 


-v' 


mixing anj 


j;le 




Fo^^(O) 




0.76 ± 0.03 










-14.5" 




The Wilson coefficients for D decays 


Clinic) 


C2{mc) 




Cl 


(1 GeV) 


C2(l GeV) 


1.216 


-0.422 


1.275 


-0.510 



annihilation corrections. To take into account the BRs and their ratio 



Ri 



(ai + a2)fn {m 



mZ 



F^- {ml 



a small l/A'J^ and Fq^'^(O) > Fq^-'^(O) are preferred (see also the discussion in footnote . 



(4.5) 



B. 



vr+TT vs. D 







K+K- 



The experiments have measured the ratio 

r(Z)° ^ K+K^) 



r(z)o 



TT+TT 



2.82 ±0.01, 



(4.6) 



which is a long-standing puzzle because the conventional factorization approach yields i?2 = 1 in 
the SU(3) limit (see discussions in Ref. [s^l). We found that the SD weak annihilation contributions 
together with FSIs interfere destructively to the — > t:'^-k~ amplitude, but constructively to the 
K^K~ amplitude, such that the ratio can be accounted for. 



In the limit of SU(3) symmetry, the — > K^K amplitude vanishes. It was explained in 
Ref. ^31] that the non-small branching ratio of this mode may be owing to long-distance FSIs. 
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TABLE II: The branching ratios in units of 10^'^: data (^Exp) 0] vs. fitted results (^fsi)- The 
individual Xi values of decay modes corresponding to the best fit are listed. For comparison, taking 
the best fit parameters and F^^{0) = 0.76 into account, we then give (i) ;BFact by means of setting 
6 = a = and neglecting the weak annihilation corrections, (ii) ^NoAnn by means of neglecting 
only the SD weak annihilation corrections, and (iii) i^NoFSi by means of setting 6 = a = 0. Note 
that — > T]'?]' is kinematically forbidden. The errors in iSpsi and xf due to the variation of 



Decay modes 



B 



Exp 



^Fact ^NoAnn B 



NoFSI 



^FSI 



xt 



38.0 ±0.9 
23.0 ±2.2 
7.7 ± 1.1 
18.8 ±2.8 



62.20 
10.39 
2.75 
2.40 



56.20 
14.05 
3.83 
3.14 



56.69 
16.69 
2.67 
16.11 



38.02 
23.83 
7.92^ 
19.82 



fO.05 
-0.09 
fO.Ol 
-0.19 
-0.17 
■0.05 
+0.17 
-0.05 



0.00 
0.14 

0.04: 

0.13 



0.04 
0.00 
+0.01 
0.06 
0.09 
0.02 
0.05 
-0.01 



D+ 
D+ 
D+ 



> K+7] 

■ K+r]' 



0.14 
0.46 
0.12 
0.11 



0.21 
0.39 
0.10 
0.12 



0.09 
0.59 
0.08 
0.19 



0.27 



+0.01 
-0.02 



0.39 ± 0.05 
0.07 ±0.00 
0.20 ±0.01 



K TT^ 



D+ 



28.2 ± 1.9 



28.15 



28.15 



28.15 



28.15: 



-0.05 
-0.03 



0.00 ±0.00 



0.138 ±0.011 



0.36 
0.03 
0.007 
0.006 



0.31 
0.06 
0.02 
0.009 



0.24 
0.03 
0.006 
0.06 



0.15 ±0.00 
0.08 ± 0.00 
0.03 ±0.01 
0.07 ±0.00 



0.68: 



-0.28 
-0.16 



D- 



2.6 ±0.7 



2.27 



2.27 



2.27 



2.27^ 



-0.03 
-0.04 



0.22 



-0.05 
-0.04 



D+ - 
D+ - 




D^- 




D^- 






K+K- 


D^- 








D^- 


A' 




rji] 




rji]' 


Z)0- 


r]'r]' 



5.9 ±0.6 

3.0 ±0.6 

5.1 ± 1.0 



11.67 
0.77 
3.47 



10.47 

2.61 

3.02 



6.77 
0.45 
4.15 



-0.29 
-0.23 



5.73: 
2.61 ±0.01 
3.31' 



-0.14 
-0.16 



0.08: 

0.42: 
3.19" 



-0.36 
-0.04 
-0.02 
-0.01 
-0.58 
-0.47 



1.38 ±0.05 
0.84 ±0.22 
3.89 ±0.14 
0.71 ±0.19 



4.73 
0.36 
4.58 


0.09 
0.09 
0.10 
0.13 





4.14 
0.71 
3.92 
0.00 
0.35 
0.15 
0.44 
0.25 




2.76 
0.22 
5.43 
0.65 
0.25 
0.01 
0.33 
1.29 




-0.03 
-0.04 
-0.01 
-0.00 
j+0.04 
-0.06 



1.37 ±0.01 

0.73^ 

3.85^ 

0.68^ 
0.68 ±0.01 

n^+°-°3 

'-'•UO_o.02 
-I -17+0.00 
-'-•-'-'-0.02 

1.95 ±0.05 




0.02 
0.26 
0.08 
0.03 



0.01 
-0.00 

0.18 
-0.11 

0.00 

0.04 
+0.23 

0.03 
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TABLE III: The Xmin/^-°-f- fitted parameters, where we obtain a twofold solution for X'^ 
which is relevant only for decay modes involving r]Q. The errors are due to the variation of F^^{0). 

Best fit results 

Xljd.o.f. i5.3tli)/5 
6 (-46 ±2)° 

a (-21 ±1)° 

Nf -21+% 
F(f''(0) 0.83 ±0.02 

K n 1 r;+0.00 

tti U.iO_Qoi 

\Xa\ 3.84 ±0.06 

arg(XA) (-138 ±3)° 



\X'a\ 


2A5to% [or 2.18 ±0.19] 




(-138 ±3)° [or (130 ±3)°] 



Nevertheless, here we conclude that D — > occurs mainly due to nonzero SD weak annihilation 

effects originating from SU(3) symmetry-breaking corrections to the distribution amplitudes of the 
kaons.^ Moreover we find of 
in Ref. 



K 



-0.15^'^Q'g'j' in the best fit, which is consistent with the result given 



331 ] but in contrast with that in Ref 



34 ] where the value is positive.^ Note that it has 



been argued in Refs. 34, 3^ that the result given in Ref. 0] is less reliable 



D. D decays involving r/ or rj' 

It should be stressed that the SD weak annihilation and SU(3) rescattering effects enter the 
amplitudes in different ways. For instance, B{D^ ^^v) arid B{D+ — > -K+rj) are mainly enhanced 
by SU(3) rescattering, whereas B{D^ — > K^-q') receives contributions mainly from the SD weak 
annihilation. This mechanism can be further tested experimentally from the relative values of the 
iT^r],iT^r]' ,r]7] and rjrj' branching ratios. 



V. DISCUSSIONS AND SUMMARY 

We have built up a simple model that the D PP decay processes go through ^^bare" am- 
plitudes followed by elastic-SU(3) rescatterings, where the bare amplitude consists of (i) the usual 
factorization amplitudes of color-allowance and color-suppression, discussed in Sec. Ill Al and (ii) 



^ In spirit, our conclusion agrees with the result in Ref. |32l ]. where the authors used the chiral perturbation 
theory to calculate the weak annihilation effects and found that the result is proportional to m^- 

3 There also exists a solution of positive af ~ 0.19, S ~ -39°, ct ~ -12°, Nf ~ -14, Xa = 2.7e*^2° ^ x'^ = 
J 7g-zi04° [oj. = 3.5e'i52°] F(f'^(0)/F(f^(0) ~ 1.05 with a larger Xmi„ ^ 11-4. 
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DO 




SU(3) rescattering-quark exchange sU(3) rescattering-quark annihlation 




SU(3) rescattering-suppressed quark exchange SU(3) rescattering-singlet exchange 

FIG. 2: Topologies relevant to D ^ Kit. The second and third rows correspond to the long- 
distance SU(3) rescattering contributions to — > K^tt'^ originating from the tree amplitude, 
where the quark exchange and singlet exchange contribute to C, the suppressed quark exchange to 
T, and the quark annihilation to the weak annihilation. The dots denote the quark fields contained 
in (y — A) fSi {V — A) four-quark operators. 



the SD weak annihilation amplitudes (H^-exchange and/or VK-annihi 
A similar model estimate was proposed by Chernyak and Zhitnisky 
amplitude followed by SU(2) rescattering for D Kir. 



ation) presented in Sec. IIIBl 
who considered the bare 



3S 



3c 



41 



In terms of quark-graph amplitudes in the diagrammatic approach |37, 
the topologies relevant to D ^ PP decays are the tree topology "T", the color-suppressed tree 
topology "C", and weak annihilation topologies (VF-exchange and/or ly-annihilation) , shown in 
the first row of Fig. [2j It has been stressed in Ref. 4]| that in the diagrammatic approach even 
though the SD weak annihilation contributions are neglected, it is still possible for that the weak 
annihilation topologies receive sizable contributions from the LD final-state rescatterings of the 
(color-suppressed) tree amplitude T (C\ as sketched in the second and third rows of Fig. |21 (some 
estimates for LD effects see Refs. 0,QI 

Moreover, it should be stressed that the SD weak 
annihilation amplitudes Aj-^^ have sizable magnitudes comparable to the factorization amplitudes 
^fac; due to the structure of (V — A) f^i {V — A) operators in the weak Hamiltonian relevant to the 
D decays, as given in (|2.16|) the SD annihilation contributions are dominated by the topologies of 
gluon emission arising from the initial state quarks of the weak vertex, whereas the contributions 
vanish in order of if the gluon is emitted from the final state quarks, i.e., the amplitudes drawn 
in Figs, nja) and (b) cancel each other. 
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We expect that the possible double counting is reasonably negligible between the LD rescat- 
terings and SD weak annihilation amplitudes due to the following three reasons: (i) The LD 
rescatterings mainly contain gluon exchanges between the two final-state mesons, as depicted in 
Fig. 121 while the gluon emission originating from the initial-state quarks of the weak vertex gives 
rise to the nonzero SD weak annihilation amplitudes, as shown in Figs.^c) and (d). (ii) The LD 
FSIs are dominated by rescatterings of the (color-suppressed) tree amplitudes which are quite dif- 
ferent from the mechanism of the SD weak annihilation amplitudes, (iii) The LD rescattering and 
SD weak annihilation contribute to amplitudes in different ways; for instance, as seen explicitly in 
Table the LS rescattering (SD weak annihilation) interfere constructively (destructively) in the 
the -D"*" — > K^7r~^ and — > vr'^vr" amplitudes. Finally, it should be noted that, in B decays, one 
may worry the double counting problem since the nonzero weak annihilation is due to the gluon 
attached to the final-state quarks in the 2{S — P) (S" -1- P) weak vertex. 

The strong phase can be generated from the radiative corrections to the weak vertex and 
the spectator interactions. Such effects were lumped into N^^ as we calculated the factorization 
amplitudes. However, since the mag nitude of l/Nf is very small obtained in our analysis, it is 
thus reasonable to neglect the resulting strong phase; choosing a real number of A'^^^, we have a 
very nice fit since Xmm/^-'^-^- ~ {^■'^-0 5)/^ negative of (or ~ 11.4/5 for positive af). In other 
words, the LD rescattering effects should be approximately absent from ^^N^^" . 

Our remaining results are briefly summarized as follows. 

• The two modest rescattering phase differences are 6 = 827 — 5s ~ —46° and o" = — 5i — 
—21°, where the a phase enters only in the ■k'^tt'-it^tt^-K^K^-K^K -Tr^rjs-rjsrjs rescattering 
subset. 

• We obtain the weak annihilation parameter \Xa\ = 3.84 ± 0.06 = 2.45^0:46 or 2.18 ± 
0.19] with a large phase (-138 ± 3)° [(-138 ± 3)° or (130 ± 3)°], where a twofold solution 
exists for X'^. 

• The K^l^ decay occurs mainly due to the short-distance weak annihilation effects, 
arising from SU(3) symmetry-breaking corrections to the distribution amplitudes of the final- 
state kaons, but receives negligible contributions from other modes via SU(3) rescattering. 

• Our results are in good agreement with the experimental measurements. The predictions for 
the branching ratios of some unmeasured modes can be used to test our model in the near 
future. 
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APPENDIX A: FACTORIZATION AMPLITUDES 



Ai^D' ^ K-7r+) = i^V^dV:,aiMml - mj,)F^'' {ml), 
AacC^'^^^V) = i9i:v^,v:,a2fK{ml-ml)F,^^m%), 
AiUD'^ ^ K%s) = i^VudV:,a2fK[cos^ml - ml)F^\mj,) 
+ sm^iml-ml,)F^'''{ml)], 

+ cos'd{ml-ml,)F^'^'{ml)], 



AuciD+^K\+) = z^VusV^a2fK{ml-ml)F,^^ml), 
AUD+^K+n'') = -i^VusV:acnfK{ml~ml)Ft{ml) 
Ai^{D+ ^ K+'ns) = i^VusV:aaifK[cosd{ml-ml)F^\rr?K) 

+ siu^{rnl-rnl,)F^'''{ml)], 

+ cos'&{ml-ml,)F^''\m],)], 



AUD^^k\+) = i-^VudV:AaiMml-ml)Fi'''{ml) 

+a2ifKiml-ml)F^^{mj,)], 

AfUD'^^K^^) = i^VusV:^a2fK{ml-ml)F^^{ml,), 
AUD^^K^m) = i^VMa2fK[cosi9iml-m^^)F^\ml) 

+ smT?(m|, - m^,)FQ'' (m|-)], 
AUD"" ^ K%) = i^VusV:aa2fK[-sm^ml - ml)F^'^{ml) 

+ cos'd{ml-ml,)F^'''{ml)], 
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+ sin'!?(m|) - mfj,)F^'''' (ml)] 
+a2[f;cos^{ml - ml)F^^ml) 

+/,nsin^(ml,-m2)Fo^-(m2,)]} 
+f^,smi»iml-ml)Fi'^ml,)], 



+ cosi?(m|, - m^j^,)FQ '' (m^)] 
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► TT^vr ) 




-> vr TT j 







-i^VudV:^a2U{ml - ml)F^^{ml), 



+/^,[sm2T?(m^ - m2,)Fo^'''(m2,) + sini?cosi9(m|5 - m2)Fo^''(m2,)]|, 



AUD^ ^ = 0, 



- /^[cosi?(m|, - m^)FQ ''{ml) + smt?(m|, - m^,)^^ (m^)] 



- /^[- sin??(m|) - m^)Fo^''(m^) + cosT?(m|, - m^,)F^'^ (ml)] 

Aac(-C'° ^ r/or?o) = iGFVudVcdO'2 

x|/,;[sm2??(m|) - m2)F^''(m2) - sin0cosi?(m^ - m2,)Fo^'''(m2)] 



(Al) 



APPENDIX B: WEAK ANNIHILATION AMPLITUDES 



Here the basic building blocks for annihilation amplitudes corresponding to Fig. ^ are denoted 
as A^^'^\P2Pi), where the superscript i{f) indicates gluon emission from the initial (final) state 
quarks, and P2 (Pi) contains a quark (antiquark) arising from the weak vertex with longitudinal 
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momentum fraction x and y, respectively, so that the building blocks read 



A{{P2Pi) = nas [ dxdy 
Jo 

A{{P2Pi) = 0. 



$P2(x)$Pi(y) 
2 

xy 



1 



+ 



y(l - xy) x^y 



A\ can be further expressed in terms of Xa as follows: 



A\ 



TTCt. 



18{Xa - 4 + ^) + 2r^r^iX^ 



3 3 
3 



-ftO 



Tra, ^18(Xa - 4 + ^) + 2r^r^2x2 

+18af (XA + 29 - 37r2)l ; 
Tra, [i8{Xa - 4 + ^) + 2r^r^iX2 
-54af(X^ + |-^) 



iiP2 = K-,K ,Pi =7r,r?8,o, 
ifP2 = 7r,r?8,o,A = i^+,i^°, 
ifP2 = i^+,i^°,Pi=7r,r?8,o, 



Tras 
7ra., 



18(Xa - 4 + ^) + 2r^V^iX2 J ; if P^, Pi = tt, r?8,o , 

18{XA-4+^)+2ri^r^^Xl 
-ISaf (XA + 29 - 37r2)j ; 



-ISaf (4Xa + 33-4^2) 

+bA{af f{XA - 71 + 
Tra, [i8(Xa - 4 + ^) + 2{r^fXl 
+18af (4XA + 33 - 47r2) 
+54(af)2(XA-71 + 77r2) 



if P2 = 7r,?78,o,i'i = i^-,i^ 



if P2 = Pi = , 



if P2 = K- {K ), Pi = K+ {K^ 



(Bl) 
(B2) 



(B3) 



where Xa — ^ -X"^ for processes containing The complete weak annihilation amplitudes are given 

by 



Ar,iK\') 
Ats(K%8) 



-i^fDfnfK^c,V^dV:,A\{K\^), 
i%fDUfrjs%iVudV:A-2A[{r^sK') + A\(k' rjs)}, 



2V3 

Ar^iWno) = i^fDUU^ciVudV:,{A\{rioK'')+A\{Wrio)} 



7O. 



7O 
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Ar,(K°7r+) = i^fj,f,f^^c2VusV:,A\i7r+K% 
Ar,{K+A = i^fDMK^C2VusV:,A\in'^K+), 

Ar^{K+r,s) = i^fDUU^C2VM{A\{V8K+)-2A{{K+r,s)}, 
Ar^{K+r,o) = i^fDUU^C2VusV^{A\{r^oK+) + A\{K+r^o)}, 

^ThI^V) = 0, 



Ar^K+TT-) = i^fj,f^fK^ciVusV^A{{7r- K+), 

Ar,{K%s) = i^fDUr,s^c,VusV:,{A{ir,sK'^)-2A{iK%s)}, 
Ar^{K%o) = i^fDUU^ciVM{A\{voK^) + A\{K%o)}, 

^r,(7r+7r°) = ^^/^/^^csK^KH^Utt" 7r+) - Al(7r+ 7r°)} = 0, 
Ar^iTT+ris) = ^^fDUU^C2VudV:a{A\{^^+r,s)+A\{v8^^+)}, 
Ar.iK+K") = i^fj,f^^c2VuMA\{K+K'), 



ArH(i^+i^-) = z^fDfK^c,VusV:,A\{K-K+), 

ArA^'m) = -i||/z././.s^ciFwl^;,{Al(7r%8) + Al(r/8 7r°)} 

ArH(Ao) = -i^/D/./,„^ciK,y;,Ki(7rO%)+AU%vrO)}, 

^TgCrysr/o) = i-^fofvafvoj^ciiVudV^ - 2VusVcs){A\{rio rjs) + ^U??8 ??o)}, 

^TbMo) = i^fDfS^c,iV^dV:i + VusV:,)A\{mVo). 
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The above weak annihilation amplitudes can be further expressed in terms of X^2 ^ follows (in 
units of i^fofpJPij^'^Oisy. 



18{Xa - 4 + y ) + 2r^rJXl + 54af (X^ + ^ - 



4 TT^ 



1 

72 



Ar^iK-7r+), 



1 



18{Xa - 4 + ^) + 2r^rfXl - 18af (5Xa + 62 - Ttt^) 



TT 



18{X'a - 4 + ^) + 2r^r'T°X'2 + 9af (2X;i - 25 + 2tt 



K I 



1 



ArsiK+A 



ArAK-^Vs) 
ArAK^m) 



18{Xa - 4 + — ) + 2r^rjXl + 18af (X^ + 29 - Stt^) 



V2 



"3 

,2 



18{Xa - 4 + ^) + 2r^r^^Xl - 18af (7Xa + 37 - bn^) 



18iX'A - 4 + ^) + 2r^rl°X'2 - 9af {2X'a - 25 + 27t^) 



ArAWTr+) = 0, 



AtAK^tt' 



1 V^sV^c, 



Ars{K%o) 



V2V^dV*C2^^'^ 

udV*sC2 



ArAK\+), 



VudV*C2 



VudV*sC2 



Ar^K+rio), 



Ar^TT+TT^) = 0, 



Ars(7r+r]s) 



C2VudV* 



18{XA-4+—) + 2rlrl^Xl 



Ar^K+lT) = C2VM 



18{Xa - 4 + y ) + 2{r^fxl 



ISaf (4Xa + 33 - 47r^) + 54(of - 71 + 7n^) 



K\2, 



Arsin^Vo) = -^c2VudV*a 



18{X'^-4+-) + 2rlrl^X'2 
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Ar^iTT+TT ) = ClVudV^ 



1 
72 



Ar^irism) = 
Ar^{K+K-) = c^VusV, 



18(X^-4+^) + 2(r^)2xl 



+18af (4X^ + 33 - 47r^) + 54(af - 71 + Ttt^) 



K\2i 



ATsi-^^m) = T^ciVudV*^ 



V3 



^rB(7r°r7o) = - 
Ar^{r]^m) = - 



18(X^-4+y) + 2(r;j)2xl 
18(X^-4+y)+2r-r^°X:f 



(B4) 

where Xa is treated as a universal parameter for SU(3) channels, while for decay modes involving 
rjo, it is distinguished to be X^. 

APPENDIX C: SU(3) FINAL STATE INTERACTIONS — 8 8 DECOMPOSI- 
TION 









fu\ 


<l = q' = 






d 











To describe elastic SU(3) final state interactions among D — > P1P2 decays, we adopt the nota- 
tions: 



(CI) 



and 

q = 9j = ( gi 92 qz) = {ud s). (C2) 
The octet final-state pseudoscalar mesons Pi and P2 , which are viewed as composites of quarks in 
the quark model, can be represented by the matrix 



n = q (g) q - -1 Tr(q q) 



V 



\/2 V6 

7r~ 



(C3) 



where W- is the 8 representation, while 11* = 0. The SU(3) final-state rescatterings for D P1P2 
are described by the product 8 8. Since the P1P2 states obey the Bose symmetry, only the 
symmetric states given by the representation 36(= 27 © 8 © 1) in 8 (X" 8(= 36 © 28) decomposition 
are relevant, whereas states given by the representation 28(= 10 © 10 © 8) vanish. The weak decay 
amplitudes Af^^ for D ^ P1P2 with FSIs are given by 

(C4) 



, FSI 



j;sf A^- = (UTSV2,,,,U),A 



bare 
I ' 
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where A}"^^^ = Ap= + aJ'^ are defined in Eqs. (EUl)) and (nTTTl) . In orthonormal 

1 /2 

bases of SU(3), the S^- matrix, describing the SU(3) FSIs, can be recast into the fohowing form 



.1/2 
diae 



27 



e*^-^ IT (27) ; a) (T (27) ; a| + e^^«^ |T (8) ; 6) (T (8) ; 5| + e^'^ \T (1)) (T (1)| , (C5) 



a=l 



where |T(27); a), |T(8); 6), and |T(1)) are orthonormal SU(3) bases in the irreducible representation 
36. Using the tensor approach 



Jll)lai 



6=1 
inorr 

the 36 states are described by n^^^II-'^^j^ with being 
symmetric in indices and can be decoupled into three types of irreducible tensors: (i) 1, an 
irreducible tensor of rank (0, 0), equals to 11*^,11*^^ = T^. (ii) 8, an irreducible tensor of rank (1, 1), is 
equivalent to TJ2 - {^/^WkT]^ = n'^H^m - (l/3)<5^n'7n*„ = Ul (iii) 27, an irreducible tensor of 
rank (2, 2), is given by T^^ +T^;- (1/5) (5^r,':^^ + 5ir- + 5^^^^^^ 

V^Y • We summarized the orthonormal states in the representation 36 together with their quantum 
numbers S and I as follows, (i) In the representation 1, the normalized state \T [1)) is 



(5 = 0,1 = 0): i=(V2| 



vr+TT ) + |7r°7r°) + jr/gr^s' 



+ V2\K+K-) + V2\K^K^)) . (C6) 



(ii) In the representation 8, the normalized states |r(8) ; 6) are 
1 



{S = 1,I=-) 



10 

T 



(VedK+TT-) - V3\K\^) - \K%s)y, 



(C7) 



(5 = 0, I = 1) 



(5 = 0, I = 0) 



3|kV0) - |K%8)), 



'l{-V6i\K\~) - V3\K-7r') + \K-rjs)y, 



vr+r/s) + 



5' 



■{V3\K+K-) - V3|kV) + 2| As)), 



(C8) 



(C9) 



— (-2|7r+7r" 
10^ ' 



2|7r°7r0) + V2|r?8%> + \K+K-) + ji^V")). (CIO) 
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(iii) In the representation 27, the normahzed states \T (27) ; a) are 
(5 = 2,7=1) 



{S = -2, I = 1) 



\k'k-), ^\k-k-)- 

1 



(S = l,/=-): |K+7r+), ^(|i^07r+) - V2|K+7rO)), 
-^(|K+7r-) + V2|KV)), li^OTT-); 
(5=-!,/=^): |:^V), 4s(|K"7r+) + ^/2|KV)), 



('5 = 1, 



^/3 
1 



(|;r'7r-) - V2|i^-7r°)), |i^-7r-); 
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(2|7r+7r-) + V2|7r°7r°) + Qv^lr/gr/g) - 6|7^+ii:-) - Q\lfK^)). (C19) 



Using the above results, one can immediately obtain the relevant U matrices and the corresponding 
SU(3) eigen-amplitudes in D decays: 
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where A^-^^q, Af-^^o, Ab-^^^, A^-^^^ and AJ'-'^^o have been defined in Eqs. (jH, (ESI), (EU, (ED 
and (|3.1ip . respectively. 
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